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Abstract 

We report on the first observation of open charm production in neutral cur- 
rent deep inelastic neutrino scattering as seen in the NuTeV detector at Fer- 
milab. The production rate is shown to be consistent with a pure gluon-Z*^ 
boson production model, and the observed level of charm production is used 
to determine the effective charm mass. As part of our analysis, we also obtain 
a new measurement for the proton-nucleon charm production cross section at 
^/s = 38.8 GeV. 

I. INTRODUCTION 

Measurements of charged current (CC) charm production in deep inelastic (DIS) neutrino 
and anti-neutrino scattering have proven to be an excellent source of information on the 
structure of the nucleon, the dynamics of heavy quark production, and the values of several 
fundamental parameters of the Standard Model (SM) of particle physics [|l]. The only 
evidence to date of neutral current (NC) charm production in z/^A^ or i?^N scattering is an 
unconfirmed observation of J/ip production [0,^. Using a new Sign Selected Quadrupole 
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Train (SSQT) beam with the high energies of the Fermilab Tevatron, NuTeV is able to 
perform a sensitive search for NC charm production through detection of events with wrong 
sign muon (WSM) final states. This occurs whenever an interaction produces a single muon 
and the muon has the opposite lepton number as the neutrino beam. The SSQT produced 
event ratios of v^/ in neutrino and a v^/ anti-neutrino mode of 0.8 x 10^^ and 4.8 x 10"^, 
respectively, making possible the WSM identification. 

In SM NC interactions, charm quarks must be produced in pairs. As is the case for 
CC charm production, considerable suppression of the v^^N u^ccX rate occurs due to the 
non-zero charm quark mass rric. This measured suppression can be used to experimentally 
determine m^. If rric is a fundamental parameter of the SM, then it should have the same 
value, up to possible quantum chromodynamics (QCD) corrections, in NC and CC u^N 
charm production, and in other physical processes such as the photo-production of charm 
and the spectrum of charmonium. 

In the fixed-flavor (FF) implementation of QCD |Q, NC charm pair production can be 
attributed completely to hard scatters between the the exchanged virtual boson and a 
gluon in the nucleon sea; this process, known as boson-gluon fusion, is illustrated in Fig. 
Other implementations of QCD |^,^ envision an intrinsic charm quark parton distribution 
function (PDF) c(x, Q^), that depends on the Bjorken scaling variable x and the absolute 
value of the squared momentum transfer Q^. In addition, some have suggested that 
non-perturbative QCD effects may produce an unusually large c (x, Q"^), particularly at high 
X. In FF QCD, c (x, Q"^) ^ over the range probed by NuTeV. The validity of this 
assumption can be tested with the data. 

II. EXPERIMENTAL APPARATUS AND BEAM 

The NuTeV (Fermilab-E815) neutrino experiment collected data during 1996-97 with 
the refurbished Lab E neutrino detector and a newly installed Sign-Selected Quadrupole 
Train(SSQT) neutrino beamline. Figure ^ illustrates the sign-selection optics employed by 
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the SSQT to pick the charge of secondary pions and kaons which determine whether 
or z/^ are predominantly produced. During NuTeV's run the primary production target 
received 1.13 x 10^^ and 1.41 x 10^^ protons-on-target in neutrino and anti-neutrino modes, 
respectively. 

The Lab E detector, described in detail elsewhere [0, consists of two major parts, a 
target calorimeter and an iron toroid spectrometer. The target calorimeter contains 690 
tons of steel sampled at 10 cm intervals by 84 3 m x 3 m scintillator counters and at 20 
cm intervals by 42 3 m x 3 m drift chambers. The toroid spectrometer consists of four 
stations of drift chambers separated by iron toroid magnets. Precision hadron and muon 
calibration beams monitored the calorimeter and spectrometer performance throughout the 
course of data taking. The calorimeter achieves a sampling-dominated hadronic energy 
resolution of <Jehad/ ^had = 2.4% © 87%/ y/En ad and an absolute scale uncertainty of 
6Ehad/Ehad = 0.5%. The spectrometer's muon energy resolution (dominated by multiple 
Coulomb scattering) is a e^/E^ = 11% and the muon momentum scale is known to 6E^/ E^ = 
1.0%. With the selection criteria used in this analysis, the muon charge mis- identification 
probability in the spectrometer is 2 x 10"^. 



III. ANALYSIS PROCEDURE 



Much of the analysis procedure follows that used for a search for flavor changing neutral 



currents using the same data set |11], and further details may be obtained from the article 



describing that analysis and in Ref. 12 



A. Introduction and Data Selection 

The analysis technique consists of comparing the visible inelasticity, yms = 
Eh ad/ {Ehad + E^), measured in the z/^ and z/^ wrong sign muon (WSM) data samples 
to a Monte Carlo (MC) simulation containing all known conventional WSM sources and a 
possible NC charm signal. The NC charm signal peaks at large values of t/ms because the 
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decay muon from the heavy flavor hadron is usually much less energetic than the hadron 
shower produced in the NC interaction. The largest background, beam impurities, is con- 
centrated at low t/yis in i/^ mode due to the characteristic (1 — y)^ behavior of interactions 
of the z/^ wrong-flavored beam background. 

Events in the WSM sample must satisfy a number of selection criteria ("cuts"). The 
fiducial volume cut requires that event vertices be reconstructed at least 25 cm-Fe (cm 
of iron) from the outer edges of the detector in the transverse directions, at least 35 cm- 
Fe downstream of the upstream face of the detector, and at least 200 cm-Fe upstream of 
the toroid. Events must contain a hadronic energy of at least 10 GeV (increased to 50 
GeV for the final NC charm fit), and exactly one track (the muon) must be found. The 
muon is required to be well-reconstructed and to pass within the understood regions of 
the toroid's magnetic field. The muon's energy must be between 10 and 150 GeV, and its 
charge must be consistent with having the opposite lepton number as the primary beam 
component. Requiring that the muon energy reconstructed in different longitudinal sections 
of the toroid agree within 25% of the value measured using the full toroid reduces charge 
mis-identification backgrounds to the 2 x 10~^ level. This latter number has been verified 
using the muon calibration beam. 

B. Source and Background Simulations 

Conventional WSM sources arise from beam impurities, right-flavor CC events where the 
charge of the muon is mis-reconstructed, CC and NC events where a vr or meson decays in 
the hadron shower, and CC charm production where the primary muon is not reconstructed 
or the charm quark is produced via a interaction. Table | lists the fractional contribution 
of each background component. The relatively large beam impurity background consists of 
contributions from hadrons (including charm) that decay before the sign-selecting dipoles 
in the SSQT, neutral kaon decays, muon decays, decay of hadrons produced by secondary 
interactions in the SSQT ("scraping"), and from decay of wrong-sign pions produced in kaon 
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decays. Table |T| summarizes the relative contributions of each beam source. For this analysis 
the beam sources can be further tuned using WSM data in z/^ mode. This procedure, which 
yields in passing a new measurement of a {pN — > ccX) at ^/s = 38.8 GeV, is detailed in the 
Appendix. 

After impurities, the next largest WSM source comes from CC production of charm 
in which the charm quark decays semi-muonically, and its decay muon is detected in the 
spectrometer. The primary lepton is either an electron, which is lost in the hadrons shower, 
or a muon which exits from or ranges out in the calorimeter. The Ue beam fraction is 
1.9(1.3)% in u {u) mode, and 22% of the CC charm events which pass WSM cuts originate 
from a z/g. 

Charged current charm production produces a broad peak at high y^j^ that must be 
handled with care. The CC charm background is simulated using a leading-order QCD charm 
production model with production, fragmentation, and charm decay parameters tuned on 



neutrino dimuon data collected by NuTeV ||T9| and a previous experiment using the same 



detector [^]. Overall normalization of the source is obtained from the measured charm-to- 
total CC cross section ratio and the single muon right-sign data sample. Simulated dimuon 
events are passed through the full GEANT simulation of the detector. Figure ^ provides 
a check of the modelling of this source through a comparison of the distribution of = 
Eh ad/ {Eh ad + E^2)^ where E^2 is the energy of the WSM in the event, between data and 
MC for dimuon events in which both muons are reconstructed by the spectrometer. This 
distribution should closely mimic the expected background to the distribution in the 
WSM sample. A comparison test between data and model yields a value of 19 for 17 
degrees of freedom. 

Finally, a NC cc event produces a WSM when the c (c)-quark decays semi-muonically 
in (z/^) mode. To compute effects of fragmentation, heavy quark decay, acceptance, and 
resolution, production is simulated with a Z°-gluon fusion model H] with charm mass pa- 
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rameterQ rric = 0.5 GeV/c^ and the GRV94-HO PDF set. No corrections for the nuclear 
environment are apphed, but possible effects are considered in the systematic error. The NC 
charm quarks are fragmented and decayed using procedures adapted from the CC charm 
modelling, and the resulting WSM events are then simulated with the full detector MC and 
processed with the data reconstruction code. 



IV. RESULTS AND INTERPRETATION 
A. Fits to Data 

Binned maximum likelihood fits are performed to the measured neutrino mode y^is dis- 
tribution using a model consisting of all conventional WSM sources described and a possible 
CC signal. The fitter varies the NC charm contribution in shape and level by allowing the 
charm mass parameter, rric, to float; it also varies the normalization of the beam impurities. 
Figure H shows the yms distribution for the data with the background plus fitted NC charm 
signal superposed. The shape indicates a preference towards including the NC charm signal, 
and the fit yields = 1.4211];^ GeV/c^ with a beam normalization of 1.00 ± 0.06, where 
the errors are purely statistical. The fitted value of the beam normalization validates the 
beam impurities model and the WSM tuning procedure described in the Appendix. 

In Fig. 1^ the yms distribution of WSM's is shown with the additional requirement that 
Ehad be larger than 50 GeV. This cut removes 85% of the beam impurities while keeping 75% 
of the NC charm events. Performing the fit again on this reduced sample with the background 
normalization fixed at 1.0 yields a consistent value for the charm mass of ttIc = lAOli^'^l 
GeV/c^, with the error again purely statistical. 



^The choice of a low input MC mass allows for consistent re- weighting of the cross section over a 
wide range in the final fit procedure. 
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B. Systematic Errors 



Estimates of the systematic uncertainty are obtained by varying parameters associated 
with background and signal simulations and event selection criteria within known bounds. 
Systematic errors are assumed to be independent and thus can be added in quadrature. More 
details on systematic error studies may be found in Refs. |TT| and |T^. The most important 
contributions are from modelling the level of CC charm events which reconstruct as WSM, 
from energy and momentum calibration uncertainties, from charmed quark fragmentation, 
from the choice of the gluon PDF used in the cc production model, from possible nuclear 
effects, and from the beam impurity model. 

The number of dimuon events which reconstruct as a WSM can be normalized either 
from the total right sign muon sample and the measured charm production fraction or 
from the observed number of events reconstructed with two events in the toroid. These 
two normalizations disagree by 3%, and switching to the latter normalization shifts by 
+0.10 GeV/c^. Replacing the drift-chamber-tracking-based method to reject events with two 
muons by a calorimeter-pulse-height-based algorithm leads to a further shift 5mc = -1-0.04 
GeV/c2. 

Shifting by ±1% calibration uncertainty changes rric by '^q'oI GeV/c^. Shifting Ehad 
by ±0.5% changes rric by ±0.01 GeV/c^. A total systematic error of 0.05 GeV/c^ is thus 
attributed to calibration. 



The fragmentation model, based on CC charm analysis of the same experiment [|l9 



IS 



tested by using events with the same production kinematics and the Lund string fragmen- 
tation model |2y]. This change increases by 0.14 GeV/c^. 

The only possible relatively-unknown input to the boson-gluon fusion cross section model 
besides the charm quark mass is the gluon PDF. Changing from the GRV94HO set to 
CTEQ4M raises rric by 0.04 GeV/c^. Varying the CTEQ gluon PDF according to the 
prescription given by its authors |2T| produced a maximum variation in rric of —0.04 GeV/c^, 
and this is used as the systematic error on the gluon PDF. 
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It is unclear whether the EMC correction p2[ for nuclear effects should be applied to 



boson-gluon fusion processes, so the NC charm fit is performed with and without it. Applying 
the EMC correction increases the measured rric by 0.12 GeV/c^, and this shift is included 
as a possible systematic error. 

The final systematic error is due to the size of the beam impurities. The beam fit 
described in the Appendix returns a normalization value and error for each of five sepa- 
rate beam sub-sources. To examine the sensitivity to each individual source, each source 
normalization is fixed one sigma high and low of its best fit value, and the other sources' 
normalizations are extracted. These alternative settings are then applied to u-mode beam 
impurities, and is re-extracted. The largest change occurs for scraping and the second 
largest for beam-produced charm. The sum, in quadrature, of all changes is a shift in rric of 
0.13 GeV/c2. 

The sum of all systematic errors in quadrature is 0.26 GeV/c^. 

V. SUMMARY AND CONCLUSIONS 

Wrong sign muon data in z/^Fe scattering data show clear evidence for NC open charm 
production. The result of a boson-gluon fusion fit with the GRV94HO gluon PDF set to 
the data yields rric = l-40l:[3|^ ± 0-26 GeV/c^. This value of charm mass corresponds to a 
production cross section a {v^N u^ccX) = {0.2lto\l) fb at an average neutrino energy 
(E) = 154 GeV. The rric governing NC neutrino charm production is consistent with the 
value obtained from CC neutrino charm production from photo-production, and from 
charmonium spectroscopy p7[|. 



Differential cross sections computed with this charm mass, the GRV94HO gluon PDF, 
and the gluon fusion model are compared to electro-production data in Fig. |. Our data 
is sensitive in a region that overlaps the EMC experiment ||28|, and it extends to slightly 
higher and slightly lower x. It is consistent with this electro-production data within rather 
large errors, providing evidence that the boson-gluon fusion process is probe- independent, 
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as expected from QCD. 

Finally, since the NC charm signal can be adequately described by the boson-gluon fusion 
diagram, there is no evidence for the existence of either a perturbative or non-perturbative 
intrinsic charm sea from NuTeV data. 
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APPENDIX A: DETERMINATION OF CHARM PRODUCTION IN PN 
SCATTERING AT 800 GEV FROM MODE WSM DATA 

Beam impurities are responsible for over 80% of the WSM's in u-mode. The beam 
impurities are due to scraping, hadrons (including charm) that decay before the sign- selecting 
dipoles in the SSQT, neutral kaon decays, muon decays and K tt ^ fi. Several of these 
sources are not well-constrained by previous measurements; and u-mode WSM's can be used 
to improve knowledge of their normalization. For the case of WSM's from charm decay, this 
effectively amounts to performing a new measurement of cr [pN — > ccX) . 

Secondary vr's and K^s are modelled by Malensek's parameterization of Atherton's 



data [|T^. The interaction of secondaries with beam elements (scraping) is modeled by 
GHEISHA [^. The production is handled by extending Malensek's charged kaon pa- 
rameterizations using the quark counting relation = {3K^ + K^) /4. The K ^ n 
process is correctly modeled. Muon decay is a well understood process; however, there 
is some uncertainty in the polarization V of the beam which should lie within the range 
V =0.1 ^Q^. Two experiments measure the inclusive cross section for production of D"^ and 
/ mesons with an 800 GeV proton beam(Ammar et al. Jl^ and Kodama et al. [0]). 



The weighted average of their production parameters are used as the starting value in this 
analysis. 

A model of the WSM's is constructed from the beam sources added, with adjustable 
weights, to non-beam sources (CC and NC charm, vr/fC decay in the shower, and charged 
mis- measurement), and a binned likelihood fit is performed jointly to the neutrino energy 
distribution and the vertical position of WSM neutrino interactions in mode . The fit 
constrains the weight of beam sources, other than charm production, to be consistent with 
1.0 within their estimated a priori errors. Table |ITT| lists these errors and gives the fit results. 
The only significant deviations from unity of any source normalizations are in the scraping 
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and charm contributions. The charm result indicates that J2i <^{Di) x BF{Di u^) = 
(9.8 ± 2.2) yub be increased of almost 50% over the a priori estimate. Recalling that there 
are two mesons for each cc pair, using BF{c ^ fi) = 9.9 ± 1.2% [Q, and assuming linear 
A dependence, one obtains a{p + N —>■ cc) = (49 ± 11) fib. Figure |^ shows the agreement 
between neutino energy (E^) of data and MC before and after the fit. 

The dominant systematic error, on a{p + N ^ cc), is 5.0 /ib due to the uncertainty in 
BF{c — > fi). The only other large systematic error is 2.3 fih due to the different methods of 
rejecting events with two muons described in the systematic errors above. Systematic errors 
due to the normalization of non-beam sources, energy calibrations, parameterization of pt 
and Xf, and the use of the i/ms < 0.5 cut are small. The total of all of these sources is 5.6 
nb yielding the final result of a{p + N ^ cc) = (49 ± 11 ± 6) /ib. 



Using PYTHIA's [|20| fragmentation of c quarks into mesons one can transform Kodama's 



and Ammar's measurements into the measurements of a{p + N ^ cc) found in Table |V. 
NuTeV's measurement is consistent with these previous measurements, and has smaller 



errors. 
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FIGURES 
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FIG. 2. Schematic of the SSQT beamline 
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FIG. 3. Comparison of data to MC of Ehad/ {Ehad + ^^2) for dimuon events where both muons 
are toroid-analyzed. 



16 



50 



40 



30 



20 



10 



data 

Before fit 
After fit 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Yvis data 

FIG. 4. Distribution of yi,is in f^-mode WSM's for data(solid), background (dashed), and 
background plus fitted NC signal(dotted) . 
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FIG. 5. Distribution of for WSM's for data(solid), backgrounds (dashed), and background 
plus NC signal(dotted) with an additional requirement Ef^^ad ^ 50 GeV. 
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FIG. 6. F2^°'^"^ as a function of x, for various Q^. The bands are the gluon-boson fusion cross 
section using rric = IAO'^q'^q and the GRV94HO gluon PDF. Data points are from charged lepton 
scattering from refs. and |2^. Our data are sensitive in a region that overlaps EMC but extends 
to shghtly higher and shghtly lower x. 
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FIG. 7. The Eu distribution of WSM events in z^-mode before and after the fit. 
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TABLES 



Source z^-mode(%) ^'-mode(%) 

Beam Impurity 67 83 

Charged Current Charm 19 8 

Charge Misidentification 5 5 

Neutral Current Charm 5 2 

Neutral Current -k/K decay 2 1 

Charged Current -k/K decay 1 1 
TABLE I. Percentage of WSM's for each source in a given mode. 

z^-mode i^-mode 

Ez/ > 20 GeV Ei^ > 20 GeV 

scraping 53% 24% 

charm 10% 25% 

KO 12% 16% 

other prompt 9% 22% 

muon decay 11% 11% 

K^TT^ H 5% 2% 



TABLE II. The percentage of beam impurities due to a given source in each mode. 
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TABLE III. Results of v-mode beam fits. * The a priori charm error is not used to constrain 
this fit. 





A/a 1 n p 


Error 




pin arrn 


1.47 


0.33 


0.30* 




1.01 


0.29 


0.20 


scrape 


1.22 


0.34 


0.40 


other 


1.00 


fixed 


0.03 


muon 


0.95 


0.11 


0.07 


prompt 


1.02 


0.21 


0.10 



TABLE IV. Previous charm meson production cross-sections transformed into charm quark 
production cross-sections. 



Exp 


a{p + p ^ cc) from measurement 


a{p + p ^ cc) from measurement 


Kodama(1991) 
Ammar(1988) 


75±18±28 
51±8±13 


47±4±16 
27l^^±7 
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